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Abstract 



A model for phase transitions initiated at nucleation sites on grain boundaries is proposed and 
I""! tested against numerical simulations: a mean field approach allows to explicitly consider the granu- 

"o 

CZ3 lar structure, yielding accurate predictions for a wide span of nucleation processes. Comparisons are 

^ made with classical models of homogeneous as well as heterogeneous nucleation. While it is shown 

. that the standard theory of homogeneous nucleation fails at predicting the nucleation kinetics from 

g the interfaces, it is additionally shown that the accuracy of the usual model of heterogeneous nu- 
cleation diminishes as the nucleation rate increases, while that of the presented model is enhanced. 
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Recrystallization is a mechanism of great scientific and technological importance, en- 
countered during the thermomechanical processing of various materials including metals 
pp. The first model that efficiently captures the main features of recrystallization, namely 
the JMAK model, states that grains nucleate from points of random locations, and that the 
grains grow at a constant velocity, until impinging other neighboring growing grains. Thanks 
to its straightforwardness, this model was also used in many other situations, provided that 
the hypotheses of random nuclei and independent growing zone are met: combustion j2], 
particle physics [3j, evolution of damage under dynamic tensile loadings, to mention only a 
few, making the JMAK model a much encountered approach. 

Random distribution of nuclei is however a rather crude hypothesis. To give examples, 
materials experience damage by crack nucleation preferably at grain boundaries [HE], com- 
bustion of solid energetic materials starts at preferred sites [6j, nucleation of grains during 
recrystallization appears at prior grains frontier [7] , and crystallization can be influenced by 
confinement in a porous media [8J or contact with grain boundaries of an other material [H]. 
Thus, the problem of nucleation from interfaces is encountered in a variety of fields [TUHH] 
and the importance of structured nucleation sites - inner nucleation free volumes bounded 
by interfaces - as well as grain size dependence is commonly witnessed [ T5\ [16]. 

In this regard, extensions of the JMAK model have been proposed over the years to take 
into account the specificities of structured nucleation. Most derivations still consider random 
distribution of nuclei, improving only marginally the model, by fitting a so-called JMAK 
parameter |17) with no clear physical justification [TBl - 121) . 

One of the major improvement of the JMAK model that faces the problem of nucle- 
ation heterogeneity has been proposed by J.W. Cahn and considers nuclei distributed on 
planar interfaces p!5| [22]. This approximation extends the predictions to nucleation at grain 
boundaries, provided that the density of nucleation sites remains low enough. Cahn's main 
assumption is that the superimposing planes are randomly located, whereas an accurate 
modeling for the higher nucleation densities requires to capture the deterministic nature of 
the location of nucleation -the grains boundaries-, excluding grain volume as a possible nuclei 
source. In other words, an assembly of random planes is a coarse description of interfaces 
in a granular material. 

In this Letter, we propose an accurate modeling of the kinetics of grain nucleation and 
growth that takes explicitly into account grain boundaries as preferred sites for nuclei. A 
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FIG. 1: Upper picture: Schematic representation of the horizon of a point g in a time-space diagram. 
Nucleation points are represented as red dots, with their associated (shaded) cone of transformation. The 
grain is represented by the two green Hnes. Lower picture: Schematic representation of the two sets of 
nucleation sites within a shce of H{q,t) at a given nucleation time t' , on interfaces (closed red dots) with 
intensity as(t') — a{t') and in volume (open dots) with intensity av{t') = y^a{t'). Point q can be 
transformed only by a nucleation inside its horizon H{q,t). 

characteristic length L„ is introduced and compared to the average grain size Lg, allowing 
for a detailed analysis of the interplay between length scales of initial microstructure and 
reacted one. It is shown that for high nucleation rates, the reacted volume is dictated by 
the granular structure only, leading to a deterministic behavior of the nucleation process. 
On the other hand, for a decreasing site density, the effect of microstructure progressively 
vanishes whereas the probabilistic nature of the nucleation process increases. The validity 
of this modeling is supported by simulations of reactions (phase transition, combustion, 
damage, ...) initiated by random nucleation on the interfaces of a Voronoi tessellation. 
Once nucleated at time t', a transformation zone expands from the nucleation site over a 
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maximum travel distance given by ht'{t) = J^* c{s)ds, c{s) being the expansion celerity and t 
the present time. All potential nucleation sites in this expansion zone then become inhibited. 
Thus, a nucleation point q' is inhibiting all nucleation in a time-growing transformed zone 
of radius hti{t). 

Considering a random point a nucleation event taking place at time t' at a distance from 
it smaller than hf (t) is close enough to be certain to transform q before or at time t. Therefore 
the probability P3 that q is transformed at time t by a nucleation that previously occurred 
is given by the probability of at least one nucleation occurring in its horizon H{q,t) = 
{{x,t'), \\x — q\\ < ht'{t)} consisting in the time-space that can influence the considered 
point q. The probability of at least one event P3 is complementary to the probability of 
no nucleation event occurring and will be determined this way. As P3 is expressed 
for any point, this probability is identified as the ratio of transformed volume at time t, 
= = P^(t). A schematic representation of the time horizon in a time-space 

diagram is proposed in the upper picture of Fig. [T| 

Considering an untransformed element of volume Aq over a short time interval At 
beginning at tj, and taking them small enough so that the nucleation rate a can be considered 
as constant, then nucleation will occur in this element {qi,tj) with a probability given by: 

P{qi,tj) = a{qi,tj)AqAt (1) 

The probability that no nucleation occurs in a horizon H{q, t) can be calculated as a 
product of probabilities over the horizon itself: 

P^{q,t)= n {l~a{q„tj)AqAt) (2) 

iqi,tj)eH{g,t) 

This expression leads to the following equation when Aq and At are infinitely small: 



$(t) = 1 - P^{q,t) = 1 - exp ( - / a{q,t')dt'dq] 

V JH{q,t) J 

= l-exp(-Ar(g,t)) (3) 

with N{q, t) built as the average number of nucleation events over the horizon H{q, t). In 
case of a transformation propagating at constant celerity c, the horizon is a set of spheres of 
radii ht/{t) = c{t — t'). If the nucleation is assumed to be homogeneous (spatially random) 
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and time independent, then the nucleation rate can be extracted out of the integral, and 
the JMAK model for homogeneous nucleation is immediately recovered, integrating over t', 
and yields for a 3D system: 



For the more general case of granular materials, the nucleation rate is known to depend 
on the local structure, and exhibit different values for nucleation occurring in the bulk or at 
the grain boundaries |22]. The following development aims at predicting the contribution of 
grain boundary nucleated phase transformation only. 

Our model is grounded on a simplified representation of an assembly of grains, composed 
of spherical grains surrounded by a homogeneous material. Inside of those spherical grain 
the probability of transformation is thus isotropic and depends only on the distance from 
the grain boundary. Considering a particular grain, nucleation can occur on its surface (with 
intensity a{t)) or beyond with an intensity approximated by the average nucleation rate over 
the whole volume, a{t)^^, as depicted in the lower picture of Fig. [l 

We will now express the fraction of transformed volume which now depends on 

spatial coordinates due to the inhomogeneity of the nucleation. Let us consider a certain 
grain of radius r and inside of it, a certain point at a distance a from the grain boundary. 
From there we can express the average number of nucleation events N{qa, t) = Nr{a, t), over 
the horizon H{qa,t), simply as a function of the distance a of Qa from the nearest grain 
boundary and of the grain's radius r. 

By noting that nucleation is not possible inside of a grain, the volume integral in equa- 
tion |3] can be split into two terms, accounting for nucleation at the surface of the current 
grain of radius r, noted 5*^ and from the rest of the material, noted The first term reduces 
to an integral over a surface defined by the intersection of the horizon with the current grain 
boundaries, and the second term corresponds to the average contribution of all other grain 
surfaces: 



In order to calculate Sr fl H{qa, t) and K H H{qa, t), we define ws, the fraction of current 
grain's surface 5*^ intersecting the horizon, which is a function of grain radius r, of the 
distance a from the grain boundary, and of the propagation time t = t — t' between a past 




(4) 
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event occurring at t' and the current time t: 



ws{r, a, t) 



if cr < a, 

TT [cV^ - a?] (^) if a < cr < 2r - a, 

47rr^ if cr > 2r — a 



(6) 



and wv the fraction of the volume K intersecting the horizon: 



if 6 < 0, 

cr + a — ISf (2cr — a + 6) 



(7) 



■p2(3r-6) ifO<6<2r, 



if 6 > 2r 



with 6 - 



2(r-^- ^ given grain radius, the investigated integral Nj.{a, t) is now expressed 
as a simple time convolution: 

t V 



Nr{a,t) 



ws{r, a,t — t') + wv{r, a,t — t 



vT 



tot 



a{t')dt' 



(8) 



The points inside of the grain and located at a distance a from the grain's boundaries de- 
fine a surface of constant probability of transformation [1 — exp {—Nr{a,t))], corresponding 
to the transformed fraction of this surface. As the probability of transformation of a point 
in a grain depends solely on its distance from the grain boundary a, we can now express 
the transformed fraction of a grain ^g{t,r) as an integration towards its center following a 
homothetic path: 

2 



[l-exp(-Nr(a,t))] 



r — a 



da 



(9) 



Thus we determined the transformed fraction of a grain ^g{t,r) given its radius r and are 
now able to express the global transformed fraction as an integration over the grain 
size distribution /(r). 



^g{t,r)f{r)dr 



(10) 



However approximated, the presented model is exact for two important limit cases. For 
very low nucleation rates a{t), the transformation of a point is most likely to be caused by 
a nucleation outside of the current grain, and N becomes significant only for time ct ':$> r, 
when ws — 47rr^ <S wy — ^[(i — t')c]^. The transformed fraction then reduces to 
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$(t) = 1 - exp{-N) with = ^ a{t')[c{t - t')]^dt', the classical solution of the JMAK 
model. 

Contrarily, a nucleation rate high enough ensures that [1 — exp(— A^)] = 1 for a < ct 
(and otherwise), yielding the exact expression of the homothetic transformations of an 
assembly of regularly shaped grains once injected in Eq 9 = ^ J^^^([r^-l)/rf da = 

[1 — [(r — min(ct, r))/r]^]. 

This model allows for a continuous representation of nucleation from heterogeneous to 
homogeneous cases. If the grain boundaries are saturated with nucleation events the mi- 
crostructure of the system dictates the kinetics of the transformation. On the contrary, if the 
nucleation events are rare enough they are not spatially correlated and appear be randomly 
distributed. We will now examine the condition of transition between these two extreme 
behaviors. 

We propose to define a characteristic time tc as the time required for one nucleation 
event to happen on average in a characteristic zone which extent is given by Z{t) = kc^t"' 
for dimension n {n = 2 for site saturated nucleation, n = 3 for diffuse nucleation), with k a 
shape parameter (/c = vr in 2D and A; = ^ in 3D). The previous condition of one nucleation 
event on average, for a constant nucleation rate: 

ankd\tc - t'Ydt' = 1 (11) 



leads to the characteristic time: 



t. 



n + l 



1 



n + l 

(12) 



ankc 

with a2 = <y and = aStot/Vtot (nucleation on surface averaged over the volume). The 
radius of the characteristic zone at tc defines a characteristic length L„, function of space 
dimension n, representing the average distance between nucleation sites: 



c(n + l] 



1 

n + l 



(13) 



kttn 

For n = 2, L2 should be orders of magnitude smaller than the grain size Lg to guarantee 
the 2D nucleation hypothesis. In the other hand, for n = 3, homogeneous nucleation can only 
hold when L3 is compatible with an averaging over numerous grains, that is for ^ Lg. A 
possible transition definition is to consider the frontier between the two domains, at L2 = L3 
leading to a definition of a transition scale Ltrans = ( ^ 
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Nucleation can thus be defined as homogeneous or heterogeneous depending on how the 
characteristic length L„ compares to this transition scale. If L„ is smaller than Ltrans? the 
characteristic zone containing one nucleation site is fixed by the granular structure and the 
heterogeneity of the system is of importance to its dynamics. On the contrary, if L„ is of the 
same order of magnitude as Ltrans, the characteristic zone containing one nucleation point is 
not imposed by the granular structure and the system can be considered as homogeneous. 
This transition is observed on Fig. |2| where the error between simulations results and 
different models are plotted over the ratio L3/L2, showing that the presented model yields 
accurate predictions over the entire domain. 



1 


1 ' 1 ' 


1 ' 1 ' 1 ' J 
/ 










■■ Simulation 


y 




- - Presented model 


/ 




■- JlVrAK 


/ 




- Calin 


/ 






/ 

/ 

y 










^ — — 

- ^ — 

1 1 1 1 1 , 1 , 1 , ] , J 




0,8 0,9 1 


1,1 1,2 1,3 1. 




L3/L2 



FIG. 2: Simulations of three-dimensional transformations nucleated from the interfaces (Voronoi tesse- 
lation) of a granular structure. Error of the models at the time at which the JMAK model predicts a 
complete transformation. The Mathematica program computing the presented model is supplied as online 
supplemental material. 

Two applications of the presented model are proposed: the kinetics of recrystallization of 
a steel, and the kinetics of crystallization of a copolymer. The relevant physical parameters 
are obtained or estimated from published results [231 - 125) : the grain size distribution / and 
nucleation rate a are derived from the analysis of the materials' microstructure, and the 
celerity of the transformation is estimated as -Rg/tioo% with Rg the average grain radius and 
tioo% the time for complete transformation. The results on Fig. |3] show that the JMAK 
model is completely unable to capture the kinetics of heterogeneous nucleation and growth 
of a transformation (unless fitted without physical justification). It also shows that Cahn's 
extension of JMAK can be very accurate in some systems (here the recrystallization of 
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steel) provided its microstructure is close to random planes. Additionally it shows that 
after inputting physical parameters, the presented model can reproduce with accuracy the 
kinetics of a transformation, demonstrating the importance of the microstructure's geometry 
as regards transformation initiated by heterogeneous nucleation. 
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FIG. 3: Experiments against models. Left: Recrystallization of an interstitial-free steel at 630°C (70% 
cold rolled IF steel in [23J). Right: Crystallization of a propylene-ethylene copolymer at 110°C (PE 5.8 in 
[21 [55]). 

We believe that the presented model has even more general applications as the mean field 
approach coupled with the space-time horizon method allows the modeling of the kinetics of 
any complex structure, Eq. [8] considering its sub-domains separately. Notably, the elongated 
form of very anisotropic grains such as those encountered in certain varieties of steel can be 
taken into account explicitly through a modification of expressions [6] and [7} Additionally, 
the kinetics of more complex nucleation processes can easily be retrieved such as a nucleation 
rate depending on temperature, through the modification of a; or type of grain, through 
the modification of a and the use of a grain distribution / considering grain variety; or a 
different rate for bulk or surface nucleation, through the addition of a bulk term with a bulk 
nucleation rate in Eq. |8j 

In conclusion, we have proposed a model for the kinetic of phase transformation nucleated 
at interfaces. We showed that the use of a mean field approach allows to reconcile heteroge- 
neous and homogeneous nucleation in a single model, and that hypotheses on the regularity 
of the granular structure make it possible to solve our problem with great accuracy even for 
very irregular structures such as Voronoi tessellations. 
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